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Abstract. - Standard models of sequential adsorption are implicitly formulated in a scale 
invariant form, by assuming adsorption on an infinite surface, with no characteristic length 
scales. In real situations, however, involving complex surfaces, intrinsic length scales may be 
relevant. We present an analytic model of continuous random sequential adsorption, in which 
the scale invariance symmetry is explicitly broken. The characteristic length is imposed by a 
set of scattered obstacles, previously adsorbed onto the surface. We show, by means of analytic 
solutions and numerical simulations, the profound effects of the symmetry breaking on both 
the jamming limit and the correlation function of the adsorbed layer. 



The irreversible adsorption of colloidal particles onto a solid surface ^] has been the 
subject of a particularly active research effort over the last years. The kinetics of adsorption 
has been mainly studied through the formulation of different models, aiming to capture the 
essential features of the process. These models are defined via a set of rules by which the 
particles accommodate when arriving at the surface, and their main purpose is to reproduce 
the experimentally observed properties of the adsorbed phase. Among those, we emphasize 
the maximum density of adsorbed particles — the jamming limit — and the structure of 
the adsorbed layer, as measured by the correlation function g(x) m. In the basic models 
proposed so far (the random sequential adsorption model (RS A) j| |[ [|, ||, 0, |§| , the ballistic 



model |9L 10, [Ll], and their subsequent extensions), the kinetics is considerably simplified, 



both numerically and analytically, by assuming that the adsorption takes place onto a planar 
surface of infinite extension. Under this condition, and given that there is no characteristic 
length involved in the problem (apart from the size of the particles) it is presumable that the 
basic dynamic quantities satisfy certain scaling laws. The existence of those scaling laws is 
ultimately responsible for the fact that, once the model is established by predicating a certain 
set of rules, all relevant quantities remain fixed and their values depend only on the rules 
adopted. In particular, all properties are independent of the particle size and the rate of 
adsorption — this latter assumed to be constant in space and time. Moreover, in the absence 
of external forces or interactions among the particles |l2|, |l3|, [l4| , the adsorbed layer exhibits 
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a simple structure in which correlations decay very fast, thus indicating absence of long range 
order. Although models formulated under those conditions may accurately reproduce many 
experimental situations, their extension to more complex surfaces having an intrinsic structure 
is by no means trivial. For example, a characteristic length may be present in the substrate, 
which may alter the kinetics of the process and the structure of the adsorbed layer. With the 
only exception of some discrete models formulated to analyze adsorption in the presence of 



point-like quenched impurities ]15| . 16 , and a model of RSA of spherical particles adsorbing 
onto a substrate composed by a random collection of points Jl7| , the physics on those complex 
substrates (beyond the scaling regime) remains essentially unexplored. 

Our purpose in this Letter is to show that when the scale invariance is broken by in- 
troducing a characteristic length, new aspects of the problem arise, leading to considerable 
changes in the adsorbed layer. The appearance of new scales may originate from the exis- 
tence of pinned objects on the surface |]l5| , |l6| . Thus, to make our analysis concrete, we will 
present an extensive study of a one-dimensional (Id) model in which particles adsorb accord- 
ing to the RSA rules onto a line where a distribution of spherical obstacles of different size 
has been previously dispersed. We note that this model is essentially different from those of 
Refs.[Ts|, 19, ^l], ^2|, which deal with the simultaneous adsorption of particles of different 



sizes. 



In the RSA model |}, the adsorbing particles are sequentially located at random po- 
sitions on the surface. If an incoming particle overlaps with a previously adsorbed one, it is 
rejected; otherwise, it becomes irreversibly adsorbed. The continuous version of this model 
can be solved in Id by analyzing the density function of gaps — holes between two consecu- 
tively adsorbed particles. In the RSA of spheres of a single diameter er, the density of gaps of 
length x, G a (x,t), fulfills the equations fic| 



I* oo 

= -{x-G)G a {x,t) + 2 G a {y,t)dy, x > a; (1) 

J x-\-a 



dG a (x,t) 
dt 

dG„(x,t) r°° 

— g- f — = J G ^(y^> d y^ x <&- (2) 

Here time t has been rescaled by the rate of arrival of the particles to the line and it has 
therefore units of inverse length. The subscript a explicitly denotes the dependence on the size 
of the particles. The coverage p a {t) is given by 1— J °° xG a (x, t)dx. The solution of Eqs. (Q) and 
(H), with the initial condition of an infinite clean substrate [i.e., G a (x, 0) = 0] yields the result 
P<r(t) = ip(ot) [|l0), where ip(t) = J Q F(u)du is the coverage corresponding to an RSA process 

with particles of size 1. We have defined the usual function F(t) = cxp |— 2 f* dz{\ — e~ z )/z j. 

From this solution we see that p^ — \im t ^oo Pa (t) = PR = 0.74759 [Q, independent of a. 
This fact can be understood by noticing that the gap distribution, as given by Eqs. (|l|) and 
(H), fulfills the identity 

G\a(x,t) = X 2 G r7 (X 1 x,Xt) 1 (3) 

for any real positive number A. This identity means that the gap density G a is scale invariant: 
covering the line (and in general any surface) with particles which are larger by a factor A 
has the only effect of reducing the time at which a given configuration is reached, by a factor 
of A -1 . It is the existence of this scale invariance that is responsible for the fact that the 
jamming limit remains fixed upon variations of the particle size a. 

The presence of impurities — defined as preadsorbed obstacles of a fixed size a® ^ a — 
breaks the scale invariance by introducing an external length scale. The jamming limit and 
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the structure of the adsorbed phase must therefore depend in this case on the size of the 
adsorbing particles, more precisely, on the size ratio of particles to obstacles, r = a/ao- 

We consider in general the adsorption of particles of size a onto a linear substrate, over 
which there is a preadsorbed set of obstacles of size ao, present with an initial density po- As- 
suming that the impurities have been adsorbed onto the surface following the RSA dynamics, 
they reach the coverage po at time to, 

po = ip(a t ) = V>(Co). (4) 

The dependence of po on to is through the combination (q = (Jota- At this time, the surface 
exhibits a gap distribution given by G ao (x, to). The problem translates then to solving Eqs. (|j) 
and (||) with the initial condition that, at time to, the gap density is given by Go{x) = 
G aa (a;, t ). 

We contemplate two possibilities, namely a > ao and a < ao- The instance a > ao 
has been previously considered in the literature, numerically |l5|, [l6| and analytically 
fib] . p3| , |24|, on one- and two-dimensional lattice (discrete) models. In this case, the imposed 
characteristic length is small, thus implying slight variations of the form of the surface with 
respect to the planar form. The first characteristic we want to show is that our model admits 
in this case a very direct solutions, which agrees with the continuum limit of the lattice models 
previously proposed. 

To show this, we perform in Eq. (jl]) the substitution G(x,t) — e~^ x ~ cr ^ t H(t). The cor- 
responding equation for H(t) is then dH/dt = (2e~ xa /t)H , which is solved with the initial 
condition H(to) — Ho- The constant Ho is determined by comparison with the initial value 
Gq(x) = Go(x > a ). We obtain 

G(x, t) = r Q t 2 e- {x - a)t F(at), x > a, (5) 

where we have defined the constant 

T _ -(a-oto)to *Wo) = .-(r-l)Co F (&) ((ft 

F(at ) ~ e F(rCo)' W 

Eq. (Q) is solved by direct integration, yielding 

G(x,t) = G (x) + T / 2ue- xu F{au)du, x < a, (7) 

J to 

where Go{x) equals Gq(x) = Gq(x < ao) or Gq(x), according to the value of x. Further 
integration of the quantity xG, in the limit t — * oo, yields the jamming limit 

Poo = Po + T [pr - tp{r(o)} ■ (8) 

This result coincides with the continuum limit obtained from the lattice model in Ref. [^3| . 

We turn now to the more complex case a < ao, in which more interesting changes are 
expected, due to the fact that the planar form of the surface is considerably altered. The 
difficulty in solving the model stems from the coupling of the solution to the initial conditions. 
We consider in particular the case Co/2 < a < ao', the case a < cto/2 can be worked out along 
the same lines. The solution of the model is found for different ranges of values of x: 

(a) x > ao'- From Eqs. (||) and (||), we readily obtain 

G (Q) (x, t) = Tot 2 e- {x - a)t F(at). (9) 
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(6) cto > x > a: From Eq. (|l]), we see that the equation is coupled to G^ a \ The solution is 
G^(x,t)=e- ( - x - a ^ t - to ^ f ° 2ue- ux F{a Q u)du + T a e-^ x -^ t f 2ue- au F(au)du. (10) 



JO Jt 

(c) a > x > <To — <t: Given G^ a \ the solution follows by direct integration: 

G (c) (x, t) = Gq (x) + T f 2ue- xu F(au)du. (11) 

(d) (To — cr > x > 0: The equation for G^ is coupled to the cases a and &. The solution is 
found to be 

G (d) (x,t) = G c (x) + 2 [ dv [ 2ue- au £ x (u + v -t Q )F{aou)du 

Jto JO 

f-t pV f-t 

+ 2T / dvC x (v) 2ue- au F{au)du + T / 2ue- [a °-' T)u F{au)du, (12) 

Jto Jto Jto 

where we have used the auxiliary function 

C x {z) = i { e - M - e-to- 0- )*} (13) 

In order to estimate the covering at jamming, we compute the integral J °° xG(x, t)dx. In 
the limit t — > oo the contribution from regions a and b vanishes, since at jamming there are no 
gaps of length larger than a particle diameter. After performing some cumbersome algebra, 
we arrive at the final expression for the total coverage, counting both particles and obstacles: 

Poo = p (l+r) + 2r [ C ° (e~ ru - 2e- u )F{u)du 
+ 2rT f 



Co 



{[1 + (1 - r)u] e~ (1 - r)u + e - ™ - 2 e - u } F{ru)du. (14) 



Equation (|14| ) determines the coverage at jamming for any value of po < PR, both explicitly 
and implicitly as a function of the initial time to- We can however determine an explicit form 
as a function of po in the limiting case £o -C 1 (small initial coverage po), by Taylor expanding 
the expression for F(t). We obtain 

Pco =Pr + (1- r)(l - p R )p + 0{pl). (15) 

That is, up to corrections of order p$, p^ grows linearly with po, with slope (1 — r)(l — 
Pr). Incidentally, we note that the same Taylor expansion in valid for the expression (0), 
corresponding to a > ctq. In this case, however, since r > 1, the jamming limit decreases 
linearly with pq. 

For larger values of po we can estimate the theoretical predictions of Eq. (|lj) by integrating 
this expression up to a very large time t. Given a value of po, the corresponding time to is 
found by numerically solving Eq. (|4[). Having obtained these two values, we perform the 
integration in Eq. (|lj). Figure |l| shows the results of such integration for different values of r, 
as a function of the initial density po- The symbols represent data obtained from Monte Carlo 
simulations of the model on a line of length L = 5000er with periodic boundary conditions, 
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Fig. 1 - Jamming limit p x as a function of the initial density po. Full lines, analytic result Eq. (|l 
filled symbols, computer simulations; dashed lines, linear approximation Eq. (|l5|). 



averaging over 100 realizations. The error bars reported are standard deviations. The dashed 
lines correspond to the results obtained in the linear approximation ()l5|). 

From Figure [l] we observe that the jamming limit is very well approximated by the linear 
expression ( |l5| ) , up to a critical density p^' , above which p^ overshoots and increases faster. 
The critical density can be seen to depend approximately linearly on the size ratio, being a 
decreasing function of r. The presence of this critical density can be understood as follows: 
The effect of a small concentration of obstacles is to essentially impose a small perturbation in 
the structure of the adsorbed phase. For r <C 1 and po <C 1, one can consider the interactions of 
particles and impurities as effectively decoupled. We are then in a situation where the particles 
saturate the surface up to a density pr, leaving free a fraction 1 — pn which is filled with 
impurities with density po. This amounts effectively to a total coverage poo pn + (1 — Pr)po, 
which indeed coincides with the exact Taylor expansion, Eq. (|l5|), in the limit r — > 0. We can 
thus interpret this case as a soft symmetry breaking, which induces at most linear corrections 
to the jamming limit. When the density po increases, the interaction between particles and 

(c) 

impurities grows larger. Beyond p , the presence of impurities radically alters the structure 
of the adsorbed phase, breaking completely the scale invariance symmetry, and inducing non- 
linear corrections to the jamming limit. 

We can assess the effects of the symmetry breaking on the structure of the adsorbed 
phase by studying the particle-particle correlation function g(x). In Figure |^ we have plotted 
g(x) as a function of the reduced length x/a, for different values of p and a fixed size ratio 
r = 0.70. Data is obtained from simulations onto a line of length L — IOOOOct, averaging over 
10 different realizations. The bin width used is 1/100. For this particular value of r, we can 
estimate from Fig. [l] pffi w 0.50. For values of po < 0.50, we observe that the correlation 
function has essentially the same shape as in standard RSA. At p = 0.50, however, we observe 
the development of a secondary peak in g(x), which eventually grows and takes over for larger 
concentrations. This peak corresponds to a majority of pairs of particles separated by exactly 
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Fig. 2 - Particle-particle correlation function for fixed r — 0.70 and different initial coverages po- 

one obstacle, such that the distance between the centers of the particles in a pair is x = a + ao 
or i/d = 1 + 1/r ~ 2.429. This value has been marked in Fig. [j] by means of a straight 
vertical line. For even larger values of the initial concentration, po > 0.70, a third peak in the 
correlations is also observed. 

In summary, we have studied the kinetics of adsorption onto a line under the conditions for 
which the density of gaps among particles is not scale invariant. The analysis of the process, 
governed by the RSA rules, has revealed the existence of new and rich phenomenology beyond 
the scaling regime, including: substantial increase in the jamming limit, which depends on 
the density of inhomogeneities introduced; the appearance of new correlations among the 
particles; and also the existence of a critical density of impurities, dividing regimes of soft and 
complete violation of the scaling symmetry To this purpose, we have proposed and solved 
analytically a continuous Id model in which particles adsorb, following the RSA rules, onto 
a surface in which spherical particles of different size have been previously scattered. When 
particularized to the case in which the imposed length scale is smaller than the size of the 
particles, our model reproduces the results obtained from a previous model considering the 
adsorption on a lattice in the presence of point-like impurities. Our results may offer new 
perspectives on what concerns modelization of the adsorption phenomena in more complex 
surfaces, as the ones having a certain roughness or an intrinsic structure manifested through 
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the existence of a characteristic correlation length. In those surfaces, found in situations of 
practical interest, the property of scale-invariance, inherent to standard adsorption models, is 
no longer satisfied. 

*** 

We thank M. A. MUNOZ for helpful discussions. The work of RPS has been supported 
by the European Network under Contract No. ERBFMRXCT980183. JMR acknowledges 
financial support by CICyT (Spain), Grant No. PB98-1258. 

REFERENCES 

[1] Bartelt M. C. and Privman V., Int. J. Mod. Phys. B, 5 (1991) 2883 

[2] Evans J., Rev. Mod. Phys., 65 (1993) 1281 

[3] Flory P. J., J. Am. Chem. Soc, 61 (1939) 1518 

[4] Renyi A., Sel. Trans. Math. Stat. Prob., 4 (1963) 203 

[5] Feder J., J. Theor. Biol, 87 (1980) 237 

[6] SCHAAF P. and Reiss H., J. Chem. Phys., 92 (1988) 4824 

[7] Senger B., Voegel J.-C, Schaaf P., Johner A., Schmidt A., and Talbot J., Phys. Rev. 
A, 44 (1991) 6926 

[8] Ramsden J. J., Phys. Rev. Lett, 71 (1993) 295 

[9] Meakin P. and Jullien R., J. Phys. (Pans), 48 (1987) 1651 

[10] Talbot J. and Ricci S. M., Phys. Rev. Lett., 68 (1992) 958 

[11] Thompson A. P. and Glandt E. D., Phys. Rev. A, 46 (1992) 4639 

[12] Pagonabarraga I., Bafaluy J., and Rubi' J. M., Phys. Rev. Lett., 75 (1995) 461 

[13] Adamczyk Z. and Warszynski P., Adv. Colloid Interface Sci., 63 (1996) 41 

[14] Pastor-Satorras R. and Rubi' J. M., Phys. Rev. Lett, 80 (1998) 5373 

[15] Milosevic D. and Svrakic N. M., J. Phys. A: Math. Gen., 26 (1993) L1061 

[16] Lee J. W., J. Phys. A: Math. Gen., 29 (33) 1996 

[17] Jin X., Wang N.-H. L., Tarjus G, and Talbot J., J. Phys. Chem., 97 (1993) 4256 

[18] Talbot J. and Schaaf P., Phys. Rev. A, 40 (1989) 422 

[19] Meakin P. and Jullien R., Phys. Rev. A, 46 (1992) 2029 

[20] Senger B., Ezzeddine R., Bafaluy F. J., Schaaf P., Cuisinier F. J. G., and Voegel 

J.-C, J. Theor. Biol., 63 (1993) 457 

[21] Adamczyk Z., Siwek B., Zembala M., and Weronsky P., J. Colloid Interface Sci., 185 

(1997) 236 

[22] Pastor-Satorras R., Phys. Rev. E, 59 (1999) 5701 

[23] Lee J. W., Phys. Rev. E, 55 (1997) 3731 

[24] Ben-Naim E. and Krapivsky P. L., J. Phys. A: Math. Gen., (1994) 3575 



